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Abstract 

The longitudinal double spin asymmetry A p x for exclusive leptoproduction of p° 
mesons, p + N^p + N + p, is studied using the COMPASS 2002 and 2003 
data. The measured reaction is incoherent exclusive p° production on polarised 
deuterons. The Q 2 and x dependence of A p x is presented in a wide kinematical 
range 3 • 10~ 3 < Q 2 < 7 (GeV/c) 2 and 5 • 10~ 5 < x < 0.05. The presented results 
are the first measurements of A p at small Q 2 (Q 2 < 0.1 (GeV/c) 2 ) and small x 
(x < 3 ■ 10~ 3 ). The asymmetry is in general compatible with zero in the whole 
kinematical range. 
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1 Introduction 

In this paper we present results on the longitudinal double spin asymmetry A p x for 
exclusive incoherent p° production in the scattering of high energy muons on nucleons. 
The experiment was carried out at CERN by the COMPASS collaboration using the 
160 GeV muon beam and the large 6 LiD polarised target. 

The studied reaction is 

H + N — fjf + p° + N', (1) 

where N is a quasi-free nucleon from the polarised deuterons. The reaction ([I]) can be 
described in terms of the virtual photoproduction process 

7* + iV -> p° + N'. (2) 

The reaction can be regarded as a fluctuation of the virtual photon into a quark- 
antiquark pair (in partonic language), or an off-shell vector meson (in Vector Meson 
Dominance model), which then scatters off the target nucleon resulting in the production 
of an on-shell vector meson. At high energies this is predominantly a diffractive process and 
plays an important role in the investigation of Pomeron exchange and its interpretation 
in terms of multiple gluon exchange. 

Most of the presently available information on the spin structure of reaction ([2]) 
stems from the p° spin density matrix elements, which are obtained from the analysis 
of angular distributions of p° production and decay [1]. Experimental results on p° spin 
density matrix elements come from various experiments [2-6] including the preliminary 
results from COMPASS [7]. 

The emerging picture of the spin structure of the considered process is the following. 
At low photon virtuality Q 2 the cross section by transverse virtual photons ot dominates, 
while the relative contribution of the cross section by longitudinal photons <jl rapidly 
increases with Q 2 . At Q 2 of about 2 (GeV/c) 2 both components become comparable and 
at a larger Q 2 the contribution of becomes dominant and continues to grow, although 
at lower rate than at low Q 2 . Approximately, the so called s-channel helicity conservation 
(SCHC) is valid, i.e. the helicity of the vector meson is the same as the helicity of the 
parent virtual photon. The data indicate that the process can be described approximately 
by the exchange in the t-channel of an object with natural parity P. Small deviations from 
SCHC are observed, also at the highest energies, whose origin is still to be understood. An 
interesting suggestion was made in Ref. [8] that at high energies the magnitudes of various 
helicity amplitudes for the reaction may shed a light on the spin-orbital momentum 
structure of the vector meson. 

A complementary information can be obtained from measurements of the double 
spin cross section asymmetry, when the information on both the beam and target polari- 
sation is used. The asymmetry is defined as 

A p x = ffl/2 ~ ff3/2 , (3) 

01/2 + 03/2 

where 01/2(3/2) stands for the cross sections of the reaction (j2J) and the subscripts denote 
the total virtual photon-nucleon angular momentum component along the virtual photon 
direction. In the following we will also use the asymmetry All which is defined for reaction 
(CQ) as the asymmetry of muon-nucleon cross sections for antiparallel and parallel beam 
and target longitudinal spin orientations. 
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In the Regge approach [9] the longitudinal double spin asymmetry A p can arise due 
to the interference of amplitudes for exchange in the t-channel of Reggeons with natural 
parity (Pomeron, p,u,f,A 2 ) with amplitudes for Reggeons with unnatural parity (-7T, Ax). 
No significant asymmetry is expected when only a non-perturbative Pomeron is exchanged 
because it has small spin-dependent couplings as found from hadron-nucleon data for cross 
sections and polarisations. 

Similarly, in the approach of Fraas [10], assuming approximate validity of SCHC, the 
spin asymmetry A p arises from the interference between parts of the helicity amplitudes 
for transverse photons corresponding to the natural and unnatural parity exchanges in 
the t channel. While a measurable asymmetry can arise even from a small contribution of 
the unnatural parity exchange, the latter may remain unmeasurable in the cross sections. 
A significant unnatural-parity contribution may indicate an exchange of certain Reggeons 
like 7r, A\ or in partonic terms an exchange of qq pairs. 

In the same reference a theoretical prediction for A p was presented, which is based 
on the description of forward exclusive p° leptoproduction and inclusive inelastic lepton- 
nucleon scattering by the off-diagonal Generalised Vector Meson Dominance (GVMD) 
model, applied to the case of polarised lepton-nucleon scattering. At the values of Bjorken 
variable x < 0.2, with additional assumptions [11], A p can be related to the A\ asymmetry 
for inclusive inelastic lepton scattering at the same x as 

< = TW (4) 

This prediction is consistent with the HERMES results for both the proton and deuteron 
targets, although with rather large errors. 

In perturbative QCD, there exists a general proof of factorisation [12] for exclu- 
sive vector meson production by longitudinal photons. It allows a decomposition of the 
full amplitude for reaction (T5]) into three components: a hard scattering amplitude for 
the exchange of quarks or gluons, a distribution amplitude for the meson and the non- 
perturbative description of the target nucleon in terms of the generalised parton distri- 
butions (GPDs), which are related to the internal structure of the nucleon. No similar 
proof of factorisation exists for transverse virtual photons, and as a consequence the in- 
terpretation of A p x in perturbative QCD is not possible at leading twist. However, a model 
including higher twist effects proposed by Martin et al. [13] describes the behaviour of 
both ol as well as of <jt reasonably well. An extension of this model by Ryskin [14] for 
the spin dependent cross sections allows to relate A{ to the spin dependent GPDs of 
gluons and quarks in the nucleon. The applicability of this model is limited to the range 
Q 2 > 4 (GeV/c) 2 . More recently another pQCD-inspired model involving GPDs has been 
proposed by Goloskokov and Kroll [15,16]. The non-leading twist asymmetry An results 
from the interference between the dominant GPD and the helicity-dependent GPD H g . 
The asymmetry is estimated to be of the order k\H,jj (Q 2 H g ), where is the transverse 
momentum of the quark and the antiquark. 

Up to now little experimental information has been available on the double spin 
asymmetries for exclusive leptoproduction of vector mesons. The first observation of a non- 
zero asymmetry A p in polarised electron-proton deep-inelastic scattering was reported by 
the HERMES experiment [11]. In the deep inelastic region (0.8 < Q 2 < 3 (GeV/c) 2 ) 
the measured asymmetry is equal to 0.23 ± 0.14 (stat) ± 0.02 (syst) [17], with little 
dependence on the kinematical variables. In contrast, for the 'quasi-real photoproduction' 
data, with (Q 2 ) = 0.13 (GeV/c) 2 , the asymmetry for the proton target is consistent with 
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zero. On the other hand the measured asymmetry A p x for the polarised deuteron target 
and the asymmetry A\ for exclusive production of meson either on polarised protons 
or deuterons are consistent with zero both in the deep inelastic and in the quasi-real 
photoproduction regions [17]. 

The HERMES result indicating a non-zero A p x for the proton target differs from the 
unpublished result of similar measurements by the SMC experiment [18] at comparable 
values of Q 2 but at about three times higher values of the photon-nucleon centre of mass 
energy W, i.e. at smaller x. The SMC measurements of A LL in several bins of Q 2 are 
consistent with zero for both proton and deuteron targets. 

2 The experimental set-up 

The experiment [19] was performed with the high intensity positive muon beam from 
the CERN M2 beam line. The fx + beam intensity is 2- 10 8 per spill of 4.8 s with a cycle time 
of 16.8 s. The average beam energy is 160 GeV and the momentum spread is u p /p = 0.05. 
The momentum of each beam muon is measured upstream of the experimental area in a 
beam momentum station consisting of several planes of scintillator strips or scintillating 
fibres with a dipole magnet in between. The precision of the momentum determination 
is typically Ap/p < 0.003. The fj, + beam is naturally polarised by the weak decays of 
the parent hadrons. The polarisation of the muon varies with its energy and the average 
polarisation is —0.76. 

The beam traverses the two cells of the polarised target, each 60 cm long, 3 cm in 
diameter and separated by 10 cm, which are placed one after the other. The target cells are 
filled with 6 LiD which is used as polarised deuteron target material and is longitudinally 
polarised by dynamic nuclear polarisation (DNP). The two cells are polarised in opposite 
directions so that data from both spin directions are recorded at the same time. The 
typical values of polarisation are about 0.50. A mixture of liquid 3 He and 4 He, used to 
refrigerate the target, and a small amount of heavier nuclei are also present in the target. 
The spin directions in the two target cells are reversed every 8 hours by rotating the 
direction of the magnetic field in the target. In this way fluxes and acceptances cancel 
in the calculation of spin asymmetries, provided that the ratio of acceptances of the two 
cells remains unchanged after the reversal. 

The COMPASS spectrometer is designed to reconstruct the scattered muons and 
the produced hadrons in wide momentum and angular ranges. It is divided in two stages 
with two dipole magnets, SMI and SM2. The first magnet, SMI, accepts charged particles 
of momenta larger than 0.4 GeV/ c, and the second one, SM2, those larger than 4 GeV/ c. 
The angular acceptance of the spectrometer is limited by the aperture of the polarised 
target magnet. For the upstream end of the target it is ±70 mrad. 

To match the expected particle flux at various locations in the spectrometer, COM- 
PASS uses various tracking detectors. Small-angle tracking is provided by stations of 
scintillating fibres, silicon detectors, micromesh gaseous chambers and gas electron mul- 
tiplier chambers. Large-angle tracking devices are multiwire proportional chambers, drift 
chambers and straw detectors. Muons are identified in large-area mini drift tubes and 
drift tubes placed downstream of hadron absorbers. Hadrons are detected by two large 
iron-scintillator sampling calorimeters installed in front of the absorbers and shielded to 
avoid electromagnetic contamination. The identification of charged particles is possible 
with a RICH detector, although in this paper we have not utilised the information from 
the RICH. 

The data recording system is activated by various triggers indicating the presence of 
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a scattered muon and/or an energy deposited by hadrons in the calorimeters. In addition 
to the inclusive trigger, in which the scattered muon is identified by coincidence signals 
in the trigger hodoscopes, several semi-inclusive triggers were used. They select events 
fulfilling the requirement to detect the scattered muon together with the energy deposited 
in the hadron calorimeters exceeding a given threshold. In 2003 the acceptance was further 
extended towards high Q 2 values by the addition of a standalone calorimetric trigger in 
which no condition is set for the scattered muon. The COMPASS trigger system allows us 
to cover a wide range of Q 2 , from quasi-real photoproduction to deep inelastic interactions. 
A more detailed description of the COMPASS apparatus can be found in Ref. [19] 

3 Event sample 

For the present analysis the whole data sample taken in 2002 and 2003 with the 
longitudinally polarised target is used. For an event to be accepted for further analysis it 
is required to originate in the target, have a reconstructed beam track, a scattered muon 
track, and only two additional tracks of oppositely charged hadrons associated to the 
primary vertex. The fluxes of beam muons passing through each target cell are equalised 
using appropriate cuts on the position and angle of the beam tracks. 

The charged pion mass hypothesis is assigned to each hadron track and the invariant 
mass of two pions, m nw , calculated. A cut on the invariant mass of two pions, 0.5 < m n7T < 
1 GeV/c 2 , is applied to select the p°. As slow recoil target particles are not detected, in 
order to select exclusive events we use the cut on the missing energy, —2.5 < E miss < 
2.5 GeV, and on the transverse momentum of p° with respect to the direction of virtual 
photon, p 2 < 0.5 (GeV/c) 2 . Here E miss = (M^-M 2 )/2M P , where M x is the missing mass 
of the unobserved recoiling system and M p is the proton mass. Coherent interactions on 
the target nuclei are removed by a cut p 2 > 0.15 (GeV/c) 2 . To avoid large corrections for 
acceptance and misidentification of events, additional cuts v > 30 GeV and Ey > 20 GeV 
are applied. 

The distributions of m nn , E miss and p\ are shown in Fig. Q] Each plot is obtained 
applying all cuts except those corresponding to the displayed variable. On the left top 
panel of Fig. [T]a clear peak of the p° resonance, centred at 770 MeV/c 2 , is visible on the 
top of the small contribution of background of the non-resonant 7r + 7r~ pairs. Also the 
skewing of the resonance peak towards smaller values of m^, due to an interference with 
the non-resonant background, is noticeable. A small bump below 0.4 GeV/c 2 is due to 
assignment of the charged pion mass to the kaons from decays of </> mesons. The mass cuts 
eliminate the non-resonant background outside of the p° peak, as well as the contribution 
of 4> mesons. 

On the right top panel of the figure the peak at E miss « is the signal of exclusive 
p° production. The width of the peak, a « 1.1 GeV, is due to the spectrometer resolution. 
Non-exclusive events, where in addition to the recoil nucleon other undetected hadrons 
are produced, appear at E miss > 0. Due to the finite resolution, however, they are not 
resolved from the exclusive peak. This background consists of two components: the double- 
diffractive events where additionally to p° an excited nucleon state is produced in the 
nucleon vertex of reaction (j2J), and events with semi-inclusive p° production, in which 
other hadrons are produced but escape detection. 

The p 2 distribution shown on the bottom panel of the figure indicates a contribution 
from coherent production on target nuclei at small p\ values. A three-exponential fit 
to this distribution was performed, which indicates also a contribution of non-exclusive 
background increasing with p 2 t . Therefore to select the sample of exclusive incoherent p° 



5 





"I Q Liiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii.il 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Pt 2 [(GeV/c) 2 ] 



Figure 1: Distributions of (top left), E miss (top right) and p 2 (bottom) for the exclu- 
sive sample. The arrows show cuts imposed on each variable to define the final sample. 



production, the aforementioned p\ cuts, indicated by arrows, were applied. 

After all selections the final sample consists of about 2.44 million events. The dis- 
tributions of Q 2 , x and W are shown in Fig. [2j The data cover a wide range in Q 2 and 
x which extends towards the small values by almost two orders of magnitude compared 
to the similar studies reported in Ref. [17]. The sharp edge of the W distribution at the 
low W values is a consequence of the cut applied on v. For this sample (W) is equal to 
10.2 GeV and (p 2 ) = 0.27(GeV/c) 2 . 



4 Extraction of asymmetry A p x 

The cross section asymmetry A LL = (a^ — cr-rt) / (^Tl + a Tt) f° r reaction ([T]) , for 
antiparallel (||) and parallel (f f ) spins of the incoming muon and the target nucleon, is 
related to the virtual-photon nucleon asymmetry A p x by 

All = D (A{ + r/A?,) , (5) 

where the factors D and rj depend on the event kinematics and A p 2 is related to the 
interference cross section for exclusive production by longitudinal and transverse virtual 
photons. As the presented results extend into the range of very small Q 2 , the exact 
formulae for the depolarisation factor D and kinematical factor rj [20] are used without 
neglecting terms proportional to the lepton mass squared m 2 . The depolarisation factor 
is given by 

n(l n 2, = y[(l + 1 2 y/2)(2-y)-2y 2 m 2 /Q 2 } 

> y 2 (1 _ 2m 2 /Q 2 ) (1 + 7 2 ) + 2(1 + #)(!-£/- 7 V/4)' W 
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Figure 2: Distributions of the kinematical variables for the final sample: Q 2 with linear 
and logarithmic vertical axis scale (top left and right panels respectively), x (bottom left), 
and the energy W (bottom right). 



where R = ctl/&Ti &l(t) is the cross section for reaction (j2J) initiated by longitudinally 
(transversely) polarised virtual photons, the fraction of the muon energy lost y = v '/E^ 
and 7 2 = Q 2 jv 2 . The kinematical factor i](y,Q 2 ) is the same as for the inclusive asym- 
metry. 

The asymmetry A p 2 obeys the positivity limit A p 2 < \^R, analogous to the one for 
the inclusive case. For Q 2 < 0.1 (GeV/c) 2 the ratio R for the reaction (j2J) is small, cf. 
Fig. [3J and the positivity limit constrains A 2 to small values. Although for larger Q 2 the 
ratio R for the process © increases with Q 2 , because of small values of r\ the product 
t]^/R is small in the whole Q 2 range of our sample. Therefore the second term in Eq. [5] 
can be neglected, so that 

K ^ ^A LL , (7) 

and the effect of this approximation is included in the systematic uncertainty of A p x . 

The number of events N{ collected from a given target cell in a given time interval 
is related to the spin-independent cross section a for reaction (T5]) and to the asymmetry 
K by 

Ni = a^ ini a(l + P B P T fDA p ), (8) 

where Pb and Pt are the beam and target polarisations, (pi is the incoming muon flux, Oj 
the acceptance for the target cell, corresponding number of target nucleons, and / the 
target dilution factor. The asymmetry is extracted from the data sets taken before and 
after a reversal of the target spin directions. The four relations of Eq. [HJ corresponding 
to the two cells (u and d) and the two spin orientations (1 and 2) lead to a second- 
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Figure 3: The ratio R = o^ja^ as a function of Q 2 measured in the E665 experiment. 
The curve is a fit to the data described in the text. 

order equation in A p x for the ratio {N u ^N d 2 /N d) iN u ^)- Here fluxes cancel out as well as 
acceptances, if the ratio of acceptances for the two cells is the same before and after 
the reversal [21]. In order to minimise the statistical error all quantities used in the 
asymmetry calculation are evaluated event by event with the weight factor w = PbJD. 
The polarisation of the beam muon, Pg, is obtained from a simulation of the beam line 
and parameterised as a function of the beam momentum. The target polarisation is not 
included in the event weight w because it may vary in time and generate false asymmetries. 
An average Pt is used for each target cell and each spin orientation. 

The ratio R, which enters the formula for D and strongly depends on Q 2 for reaction 
([2]), was calculated on an event-by-event basis using the parameterisation 



with a = 0.66 ± 0.05, and a\ = 0.61 ± 0.09. The parameterisation was obtained by the 
Fermilab E665 experiment from a fit to their R measurements for exclusive p° muopro- 
duction on protons [3]. These are shown in Fig. [3] together with the fitted Q 2 - dependence. 
The preliminary COMPASS results on R for the incoherent exclusive p° production on 
the nucleon [7] , which cover a broader kinematic region in Q 2 , agree reasonably well with 
this parameterisation. The uncertainty of Qq and cti is included in the systematic error of 



The dilution factor / gives the fraction of events of reaction ([2]) originating from 
nucleons in polarised deuterons inside the target material. It is calculated event-by-event 
using the formula 



Here % and n\ denote numbers of nucleons in deuteron and nucleus of atomic mass A 
in the target, and 5"d and a a are the cross sections per nucleon for reaction ([2]) occurring 
on the deuteron and on the nucleus of atomic mass A, respectively. The sum runs over all 
nuclei present in the COMPASS target. The factor C\ takes into account that there are 
two polarised deuterons in the 6 LiD molecule, as the 6 Li nucleus is in a first approximation 
composed of a deuteron and an a particle. 

The measurements of the <ta/5"d for incoherent exclusive p° production come from 
the NMC [2], E665 [22] and early experiments on p° photoproduction [23]. They were 



R{Q 2 ) = a (Q 2 ) 



(9) 
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Figure 4: (Left) Parameter a of Eq. [TT] as a function of Q 2 (from Ref. [24]). The exper- 
imental points and the fitted curve are shown. See text for details. (Right) The dilution 
factor / function of Q . 



fitted in Ref. [24] with the formula: 

~a A = a p ■ A a ^~\ with a{Q 2 ) - 1 = ~ exp{-Q 2 /Q 2 }, (11) 

where er p is the cross section for reaction (j2J) on the free proton. The value of the fitted 
parameter Ql is equal to 9 ± 3 (GeV/c) 2 . The measured values of the parameter a and 
the fitted curve a(Q 2 ) are shown on the left panel of Fig. H] taken from Ref. [24]. On the 
right panel of the figure the average value of / is plotted for the various Q 2 bins used in 
the present analysis. The values of / are equal to about 0.36 in most of the Q 2 range, 
rising to about 0.38 at the highest Q 2 . 

The radiative corrections (RC) have been neglected in the present analysis, in par- 
ticular in the calculation of /, because they are expected to be small for reaction (JT]). 
They were evaluated [25] to be of the order of 6% for the NMC exclusive p° production 
analysis. The small values of RC are mainly due to the requirement of event exclusivity 
via cuts on E miss and p 2 , which largely suppress the dominant external photon radiation. 
The internal (infrared and virtual) RC were estimated in Ref. [25] to be of the order of 
2%. 

5 Systematic errors 

The main systematic uncertainty of A p x comes from an estimate of possible false 
asymmetries. In order to improve the accuracy of this estimate, in addition to the standard 
sample of incoherent events, a second sample was selected by changing the p 2 cuts to 

< p 2 < 0.5 (GeV/c) 2 , (12) 

and keeping all the remaining selections and cuts the same as for the 'incoherent sample'. 
In the following it will be referred to as the 'extended p\ sample'. Such an extension of the 
p 2 range allows one to obtain a sample which is about five times larger than the incoherent 
sample. However, in addition to incoherent events such a sample contains a large fraction of 
events originating from coherent p° production. Therefore, for the estimate of the dilution 
factor / a different nuclear dependence of the exclusive cross section was used, applicable 
for the sum of coherent and incoherent cross sections [2]. The physics asymmetries A'[ for 
both samples are consistent within statistical errors. 
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Possible, false experimental asymmetries were searched for by modifying the se- 
lection of data sets used for the asymmetry calculation. The grouping of the data into 
configurations with opposite target-polarisation was varied from large samples, covering 
at most two weeks of data taking, into about 100 small samples, taken in time intervals of 
the order of 16 hours. A statistical test was performed on the distributions of asymmetries 
obtained from these small samples. In each of the Q 2 and x bins the dispersion of the 
values of A p around their mean agrees with the statistical error. Time-dependent effects 
which would lead to a broadening of these distributions were thus not observed. Allowing 
the dispersion of A p to vary within its two standard deviations we obtain for each bin an 
upper bound for the systematic error arising from time- dependent effects 

0"falseA,tdep < 0.56 fX sta t- (13) 

Here <T sta t is the statistical error on A p for the extended p 2 sample. The uncertainty on the 
estimates of possible false asymmetries due to the time- dependent effects is the dominant 
contribution to the total systematic error in most of the kinematical region. 

Asymmetries for configurations where spin effects cancel out were calculated to 
check the cancellation of effects due to fluxes and acceptances. They were found compatible 
with zero within the statistical errors. Asymmetries obtained with different settings of the 
microwave (MW) frequency, used for DNP, were compared in order to test possible effects 
related to the orientation of the target magnetic field. The results for the extended p 2 
sample tend to show that there is a small difference between asymmetries for the two 
MW configurations. However, because the numbers of events of the data samples taken 
with each MW setting are approximately balanced, the effect of this difference on A p is 
negligible for the total sample. 

The systematic error on A p also contains an overall scale uncertainty of 6.5% due to 
uncertainties on Pb and Pr- The uncertainty of the parameterisation of R{Q 2 ) affects the 
depolarisation factor D. The uncertainty of the dilution factor / is mostly due to uncer- 
tainty of the parameter a(Q 2 ) which takes into account nuclear effects in the incoherent 
p° production. The neglect of the A p 2 term mainly affects the highest bins of Q 2 and x. 

Another source of systematic errors is due to the contribution of the non-exclusive 
background to our sample. This background originates from two sources. First one is 
due to the production of p° accompanied by the dissociation of the target nucleon, the 
second one is the production of p° in inclusive scattering. In order to evaluate the amount 
of background in the sample of exclusive events it is necessary to determine the E miss 
dependence for the non-exclusive background in the region under the exclusive peak (cf. 
Fig. H]). For this purpose complete Monte Carlo simulations of the experiment were used, 
with events generated by either the PYTHIA 6.2 or LEPTO 6.5.1 generators. Events 
generated with LEPTO come only from deep inelastic scattering and cover the range of 
Q 2 > 0.5 (GeV/c) 2 . Those generated with PYTHIA cover the whole kinematical range 
of the experiment and include exclusive production of vector mesons and processes with 
diffractive excitation of the target nucleon or the vector meson, in addition to inelastic 
production. 

The generated MC events were reconstructed and selected for the analysis using the 
same procedure as for the data. In each bin of Q 2 the E m i ss distribution for the MC was 
normalised to the corresponding one for the data in the range of large E miss > 7.5 GeV. 
Then the normalised MC distribution was used to estimate the number of background 
events under the exclusive peak in the data. The fraction of background events in the 
sample of incoherent exclusive p° production was estimated to be about 0.12±0.06 in most 
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of the kinematical range, except in the largest Q 2 region, where it is about 0.24±0.12. The 
large uncertainties of these fractions reflect the differences between estimates from LEPTO 
and PYTHIA in the region where they overlap. In the case of PYTHIA the uncertainties 
on the cross sections for diffractive photo- and electroproduction of vector mesons also 
contribute. For events generated with PYTHIA the E m i ss distributions for various physics 
processes could be studied separately. It was found that events of p° production with an 
excitation of the target nucleon into N* resonances of small mass, M < 2 GeV/c 2 , cannot 
be resolved from the exclusive peak and therefore were not included in the estimates of 
number of background events. 

An estimate of the asymmetry A p for the background was obtained using a non- 
exclusive sample, which was selected with the standard cuts used in this analysis, except 
the cut on E miss which was modified to E miss > 2.5 GeV. In different h.igh-E miss bins A p 
for this sample was found compatible with zero. 

Because no indication of a non-zero A p for the background was found, and also due 
to a large uncertainty of the estimated amount of background in the exclusive sample, 
no background corrections were made. Instead, the effect of background was treated as a 
source of systematic error. Its contribution to the total systematic error was not significant 
in most of the kinematical range, except for the highest Q 2 and x. 

The total systematic error on A p was obtained as a quadratic sum of the errors 
from all discussed sources. Its values for each Q 2 and x bin are given in Tables [T] and [21 
The total systematic error amounts to about 40% of the statistical error for most of the 
kinematical range. Both errors become comparable in the highest bin of Q 2 . 

6 Results 

The COMPASS results on A p are shown as a function of Q 2 and x in Fig. [5] and 
listed in Tables Q] and EJ The statistical errors are represented by vertical bars and the 
total systematic errors by shaded bands. 




_-|L i i i i ■ ■ 1 1 ■ i ' I -|r , , ; 

icr 3 icr 2 icr 1 1 10 NT" io 3 10 2 icr 1 



Q 2 [(GeV/c) 2 ] x 

Figure 5: A p as a function of Q 2 (left) and x (right) from the present analysis. Error bars 
correspond to statistical errors, while bands at the bottom represent the systematical 
errors. 

The wide range in Q 2 covers four orders of magnitude from 3 ■ 10~ 3 to 7 (GeV/c) 2 . 
The domain in x which is strongly correlated with Q 2 , varies from 5 • 10~ 5 to about 
0.05 (see Tables for more details). For the whole kinematical range the A p asymmetry 
measured by COMPASS is consistent with zero. As discussed in the introduction, this 
indicates that the role of unnatural parity exchanges, like it- or Ai-Reggeon exchange, is 
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Table 1: Asymmetry A p as a function of Q 2 . Both the statistical errors (first) and the 
total systematic errors (second) are listed. 



Q range 


(Q ) [(GeV/c) 2 ] 


(x) 


(u) [GeV] 




0.0004 - 0.005 


0.0031 


4.0 • 10~ 5 


42.8 


-0.030 ± 0.045 ± 0.014 


0.005 - 0.010 


0.0074 


8.4 • 10" 5 


49.9 


0.048 ± 0.038 ± 0.013 


0.010 - 0.025 


0.017 


1.8 • 10" 4 


55.6 


0.063 ±0.026 ±0.014 


0.025 - 0.050 


0.036 


3.7 • 10" 4 


59.9 


-0.035 ± 0.027 ± 0.009 


0.05-0.10 


0.072 


7.1 • 10" 4 


62.0 


-0.010 ±0.028 ±0.008 


0.10-0.25 


0.16 


0.0016 


62.3 


-0.019 ±0.029 ±0.009 


0.25-0.50 


0.35 


0.0036 


60.3 


0.016 ±0.045 ±0.014 


0.5 — 1 


0.69 


0.0074 


58.6 


0.141 ± 0.069 ± 0.030 


1-4 


1.7 


0.018 


59.7 


0.000 ± 0.098 ± 0.035 


4-50 


6.8 


0.075 


55.9 


-0.85 ±0.50 ±0.39 



small in that kinematical domain, which is to be expected if diffraction is the dominant 
process for reaction (J2]). 

In Fig. E] the COMPASS results are compared to the HERMES results on A p x ob- 
tained on a deuteron target [17]. Note that the lowest Q 2 and x HERMES points, re- 
ferred to as 'quasi-photoproduction', come from measurements where the kinematics of 
the small-angle scattered electron was not measured but estimated from a MC simulation. 
This is in contrast to COMPASS, where scattered muon kinematics is measured even at 
the smallest Q 2 . 
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Figure 6: A p as a function of Q 2 (left) and x (right) from the present analysis (circles) 
compared to HERMES results on the deuteron target (triangles). For the COMPASS 
results inner bars represent statistical errors, while the outer bars correspond to the total 
error. For the HERMES results vertical bars represent the quadratic sum of statistical 
and systematic errors. The curve represents the prediction explained in the text. 



The results from both experiments are consistent within errors. The kinematical 
range covered by the present analysis extends further towards small values of x and Q 2 
by almost two orders of magnitude. In each of the two experiments A p x is measured at 
different average 14 7 , which is equal to about 10 GeV for COMPASS and 5 GeV for 
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Table 2: Asymmetry A p as a function of x. Both the statistical errors (first) and the total 
systematic errors (second) are listed. 



x range 


(x) 


(Q > l(GeV/c) z ] 


(v) [GeV] 


A l 


8 • 1CT 6 - 1 • 10~ 4 


5.8 • 10" 5 


0.0058 


51.7 


0.035 ±0.026 ±0.011 


1 • 10 -4 - 2.5 • 10 -4 


1.7- 10 -4 


0.019 


59.7 


0.036 ± 0.024 ± 0.010 


2.5 • 1CT 4 - 5 • 1CT 4 


3.6 ■ 10" 4 


0.041 


61.3 


-0.039 ±0.027 ±0.012 


5 • 10" 4 - 0.001 


7.1 • 10 -4 


0.082 


60.8 


-0.010 ±0.030 ±0.010 


0.001 - 0.002 


0.0014 


0.16 


58.6 


-0.005 ± 0.036 ± 0.013 


0.002 - 0.004 


0.0028 


0.29 


54.8 


0.032 ±0.050 ±0.019 


0.0U4 — 0.01 


O.OUoz 


0.59 


50. 1 


0.019 ± 0.0b9 ± 0.026 


0.01 - 0.025 


0.015 


1.3 


47.5 


-0.03 ±0.14 ±0.06 


0.025 - 0.8 


0.049 


3.9 


43.8 


-0.27 ±0.38 ±0.19 



HERMES. Thus, no significant W dependence is observed for A p on an isoscalar nucleon 
target. 

The x dependence of the measured A p is compared in Fig. [U] to the prediction 
given by Eq. HI which relates A p to the asymmetry A\ for the inclusive inelastic lepton- 
nucleon scattering. To produce the curve the inclusive asymmetry A\ was parameterised 
as At(x) = (x a - 7 a ) • (1 - e- px ) , where a = 1.158 ± 0.024, (3 = 125.1 ± 115.7 and 
7 = 0.0180 ± 0.0038. The values of the parameters have been obtained from a fit of 
Ai(x) to the world data from polarised deuteron targets [26-31] including COMPASS 
measurements at very low Q 2 and x [32]. Within the present accuracy the results on A p 
are consistent with this prediction. 

In the highest Q 2 bin, (Q 2 ) = 6.8 (GeV/c) 2 , in the kinematical domain of applica- 
bility of pQCD-inspired models which relate the asymmetry to the spin-dependent GPDs 
for gluons and quarks (cf. Introduction), one can observe a hint of a possible nonzero 
asymmetry, although with a large error. It should be noted that in Ref. [18] a nega- 
tive value of An different from zero by about 2 standard deviations was reported at 
(Q 2 ) = 7.7 (GeV/c) 2 . At COMPASS, including the data taken with the longitudinally 
polarised deuteron target in 2004 and 2006 will result in an increase of statistics by a 
factor of about three compared to the present paper, and thus may help to clarify the 
issue. 

For the whole Q 2 range future COMPASS data, to be taken with the polarised 
proton target, would be very valuable for checking if the role of the flavour-blind exchanges 
is indeed dominant, as expected for the Pomeron- mediated process. 

7 Summary 

The longitudinal double spin asymmetry A p for the diffractive muoproduction of p° 
meson, [i + N—^fi + N + p, has been measured by scattering longitudinally polarised 
muons off longitudinally polarised deuterons from the 6 LiD target and selecting incoherent 
exclusive p° production. The presented results for the COMPASS 2002 and 2003 data cover 
a range of energy W from about 7 to 15 GeV. 

The Q 2 and x dependence of A p x is presented in a wide kinematical range 3 • 10~ 3 < 
Q 2 <7 (GeV/c) 2 and 5 • 10~ 5 < x < 0.05. These results extend the range in Q 2 and x by 
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two orders of magnitude down with respect to the existing data from HERMES. 

The asymmetry A p x is compatible with zero in the whole x and Q 2 range. This may 
indicate that the role of unnatural parity exchanges like tt- or Aj-Reggeon exchange is 
small in that kinematical domain. 

The x dependence of measured A p is consistent with the prediction of Ref. [11] which 
relates A p to the asymmetry A\ for the inclusive inelastic lepton-nucleon scattering. 
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